The metabolic domestication syndrome of budding yeast.

Cellular metabolism evolves through changes in the structure and quantitative states of metabolic networks. Here, we explore the evolutionary dynamics of metabolic states by focusing on the collection of metabolite levels, the metabolome, which captures key aspects of cellular physiology. Using a phylogenetic framework, we profiled metabolites in 27 populations of nine budding yeast species, providing a graduated view of metabolic variation across multiple evolutionary time scales. Metabolite levels evolve more rapidly and independently of changes in the metabolic network's structure, providing complementary information to enzyme repertoire. Although metabolome variation accumulates mainly gradually over time, it is profoundly affected by domestication. We found pervasive signatures of convergent evolution in the metabolomes of independently domesticated clades of Saccharomyces cerevisiae. Such recurring metabolite differences between wild and domesticated populations affect a substantial part of the metabolome, including rewiring of the TCA cycle and several amino acids that influence aroma production, likely reflecting adaptation to human niches. Overall, our work reveals previously unrecognized diversity in central metabolism and the pervasive influence of human-driven selection on metabolite levels in yeasts.

Snowball: a novel gene family required for developmental patterning of fruiting bodies of mushroom-forming fungi (Agaricomycetes).

The morphogenesis of sexual fruiting bodies of fungi is a complex process determined by a genetically encoded program. Fruiting bodies reached the highest complexity levels in the Agaricomycetes; yet, the underlying genetics is currently poorly known. In this work, we functionally characterized a highly conserved gene termed snb1, whose expression level increases rapidly during fruiting body initiation. According to phylogenetic analyses, orthologs of snb1 are present in almost all agaricomycetes and may represent a novel conserved gene family that plays a substantial role in fruiting body development. We disrupted snb1 using CRISPR/Cas9 in the agaricomycete model organism Coprinopsis cinerea. snb1 deletion mutants formed unique, snowball-shaped, rudimentary fruiting bodies that could not differentiate caps, stipes, and lamellae. We took advantage of this phenotype to study fruiting body differentiation using RNA-Seq analyses. This revealed differentially regulated genes and gene families that, based on wild-type RNA-Seq data, were upregulated early during development and showed tissue-specific expression, suggesting a potential role in differentiation. Taken together, the novel gene family of snb1 and the differentially expressed genes in the snb1 mutants provide valuable insights into the complex mechanisms underlying developmental patterning in the Agaricomycetes. IMPORTANCE: Fruiting bodies of mushroom-forming fungi (Agaricomycetes) are complex multicellular structures, with a spatially and temporally integrated developmental program that is, however, currently poorly known. In this study, we present a novel, conserved gene family, Snowball (snb), termed after the unique, differentiation-less fruiting body morphology of snb1 knockout strains in the model mushroom Coprinopsis cinerea. snb is a gene of unknown function that is highly conserved among agaricomycetes and encodes a protein of unknown function. A comparative transcriptomic analysis of the early developmental stages of differentiated wild-type and non-differentiated mutant fruiting bodies revealed conserved differentially expressed genes which may be related to tissue differentiation and developmental patterning fruiting body development.

ContScout: sensitive detection and removal of contamination from annotated genomes.

Contamination of genomes is an increasingly recognized problem affecting several downstream applications, from comparative evolutionary genomics to metagenomics. Here we introduce ContScout, a precise tool for eliminating foreign sequences from annotated genomes. It achieves high specificity and sensitivity on synthetic benchmark data even when the contaminant is a closely related species, outperforms competing tools, and can distinguish horizontal gene transfer from contamination. A screen of 844 eukaryotic genomes for contamination identified bacteria as the most common source, followed by fungi and plants. Furthermore, we show that contaminants in ancestral genome reconstructions lead to erroneous early origins of genes and inflate gene loss rates, leading to a false notion of complex ancestral genomes. Taken together, we offer here a tool for sensitive removal of foreign proteins, identify and remove contaminants from diverse eukaryotic genomes and evaluate their impact on phylogenomic analyses.

Vertical and horizontal gene transfer shaped plant colonization and biomass degradation in the fungal genus Armillaria.

The fungal genus Armillaria contains necrotrophic pathogens and some of the largest terrestrial organisms that cause tremendous losses in diverse ecosystems, yet how they evolved pathogenicity in a clade of dominantly non-pathogenic wood degraders remains elusive. Here we show that Armillaria species, in addition to gene duplications and de novo gene origins, acquired at least 1,025 genes via 124 horizontal gene transfer events, primarily from Ascomycota. Horizontal gene transfer might have affected plant biomass degrading and virulence abilities of Armillaria, and provides an explanation for their unusual, soft rot-like wood decay strategy. Combined multi-species expression data revealed extensive regulation of horizontally acquired and wood-decay related genes, putative virulence factors and two novel conserved pathogenicity-induced small secreted proteins, which induced necrosis in planta. Overall, this study details how evolution knitted together horizontally and vertically inherited genes in complex adaptive traits of plant biomass degradation and pathogenicity in important fungal pathogens.

Genomes of fungi and relatives reveal delayed loss of ancestral gene families and evolution of key fungal traits.

Fungi are ecologically important heterotrophs that have radiated into most niches on Earth and fulfil key ecological services. Despite intense interest in their origins, major genomic trends of their evolutionary route from a unicellular opisthokont ancestor to derived multicellular fungi remain poorly known. Here we provide a highly resolved genome-wide catalogue of gene family changes across fungal evolution inferred from the genomes of 123 fungi and relatives. We show that a dominant trend in early fungal evolution has been the gradual shedding of protist genes and the punctuated emergence of innovation by two main gene duplication events. We find that the gene content of non-Dikarya fungi resembles that of unicellular opisthokonts in many respects, owing to the conservation of protist genes in their genomes. The most rapidly duplicating gene groups included extracellular proteins and transcription factors, as well as ones linked to the coordination of nutrient uptake with growth, highlighting the transition to a sessile osmotrophic feeding strategy and subsequent lifestyle evolution as important elements of early fungal history. These results suggest that the genomes of pre-fungal ancestors evolved into the typical filamentous fungal genome by a combination of gradual gene loss, turnover and several large duplication events rather than by abrupt changes. Consequently, the taxonomically defined Fungi represents a genomically non-uniform assemblage of species.

Dual RNA-Seq Profiling Unveils Mycoparasitic Activities of Trichoderma atroviride against Haploid Armillaria ostoyae in Antagonistic Interaction Assays.

Armillaria ostoyae, a species among the destructive forest pathogens from the genus Armillaria, causes root rot disease on woody plants worldwide. Efficient control measures to limit the growth and impact of this severe underground pathogen are under investigation. In a previous study, a new soilborne fungal isolate, Trichoderma atroviride SZMC 24276 (TA), exhibited high antagonistic efficacy, which suggested that it could be utilized as a biocontrol agent. The dual culture assay results indicated that the haploid A. ostoyae-derivative SZMC 23085 (AO) (C18/9) is highly susceptible to the mycelial invasion of TA. In the present study, we analyzed the transcriptome of AO and that of TA in in vitro dual culture assays to test the molecular arsenal of Trichoderma antagonism and the defense mechanisms of Armillaria. We conducted time-course analysis and functional annotation and analyzed enriched pathways and differentially expressed genes including biocontrol-related candidate genes from TA and defense-related candidate genes from AO. The results indicated that TA deployed several biocontrol mechanisms when confronted with AO. In response, AO initiated multiple defense mechanisms to protect against the fungal attack. To our knowledge, the present study offers the first transcriptome analysis of a biocontrol fungus attacking AO. Overall, this study provides insights that aid the further exploration of plant pathogen-biocontrol agent interaction mechanisms. IMPORTANCE Armillaria species can survive for decades in the soil on dead woody debris, develop rapidly under favorable conditions, and harmfully infect newly planted forests. Our previous study found Trichoderma atroviride to be highly effective in controlling Armillaria growth; therefore, our current work explored the molecular mechanisms that might play a key role in Trichoderma-Armillaria interactions. Direct confrontation assays combined with time course-based dual transcriptome analysis provided a reliable system for uncovering the interactive molecular dynamics between the fungal plant pathogen and its mycoparasitic partner. Furthermore, using a haploid Armillaria isolate allowed us to survey the deadly prey-invading activities of the mycoparasite and the ultimate defensive strategies of its prey. Our current study provides detailed insights into the essential genes and mechanisms involved in Armillaria defense against Trichoderma and the genes potentially involved in the efficiency of Trichoderma to control Armillaria. In addition, using a sensitive haploid Armillaria strain (C18/9), with its complete genome data already available, also offers the opportunity to test possible variable molecular responses of Armillaria ostoyae toward diverse Trichoderma isolates with various biocontrol abilities. Initial molecular tests of the dual interactions may soon help to develop a targeted biocontrol intervention with mycoparasites against plant pathogens.

Prolonged activity of the transposase helper may raise safety concerns during DNA transposon-based gene therapy.

DNA transposon-based gene delivery vectors represent a promising new branch of randomly integrating vector development for gene therapy. For the side-by-side evaluation of the piggyBac and Sleeping Beauty systems-the only DNA transposons currently employed in clinical trials-during therapeutic intervention, we treated the mouse model of tyrosinemia type I with liver-targeted gene delivery using both transposon vectors. For genome-wide mapping of transposon insertion sites we developed a new next-generation sequencing procedure called streptavidin-based enrichment sequencing, which allowed us to identify approximately one million integration sites for both systems. We revealed that a high proportion of piggyBac integrations are clustered in hot regions and found that they are frequently recurring at the same genomic positions among treated animals, indicating that the genome-wide distribution of Sleeping Beauty-generated integrations is closer to random. We also revealed that the piggyBac transposase protein exhibits prolonged activity, which predicts the risk of oncogenesis by generating chromosomal double-strand breaks. Safety concerns associated with prolonged transpositional activity draw attention to the importance of squeezing the active state of the transposase enzymes into a narrower time window.

Radiation of mushroom-forming fungi correlates with novel modes of protecting sexual fruiting bodies.

The protection of vulnerable developing structures evolved repeatedly in terrestrial organisms and includes, among others, viviparity in animals and the seed in land plants. In mushroom-forming fungi (Agaricomycetes), sexual spores are born on fruiting bodies, the growth of which is a complex developmental process that is exposed to environmental factors (e.g., desiccation, fungivorous animals). Mushroom-forming fungi evolved a series of innovations in fruiting body protection, however, how these emerged is obscure, leaving the evolutionary principles of fruiting body development poorly known. Here, we show that developmental innovations that lead to the spore-producing surface (hymenophore) being enclosed in a protected environment display asymmetry in their evolution and are associated with increased diversification rates. 'Enclosed' development evolved convergently and became a dominant developmental type in several clades of mushrooms. This probably mirrors spore production benefits for species with protected fruiting body initials, by better coping with environmental factors. Our observations highlight new morphological traits associated with mushroom diversification that parallel the evolution of protection strategies in other organisms, such as viviparity or the seed in animals or plants, respectively, but in the context of spore development, highlighting the general importance of protecting vulnerable progeny across the tree of life.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom-forming fungi (Agaricomycetes).

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a 'developmental hourglass,' act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Reappraisal of the Genus Exsudoporus (Boletaceae) Worldwide Based on Multi-Gene Phylogeny, Morphology and Biogeography, and Insights on Amoenoboletus.

The boletoid genera Butyriboletus and Exsudoporus have recently been suggested by some researchers to constitute a single genus, and Exsudoporus was merged into Butyriboletus as a later synonym. However, no convincing arguments have yet provided significant evidence for this congeneric placement. In this study, we analyze material from Exsudoporus species and closely related taxa to assess taxonomic and phylogenetic boundaries between these genera and to clarify species delimitation within Exsudoporus. Outcomes from a multilocus phylogenetic analysis (ITS, nrLSU, tef1-α and rpb2) clearly resolve Exsudoporus as a monophyletic, homogenous and independent genus that is sister to Butyriboletus. An accurate morphological description, comprehensive sampling, type studies, line drawings and a historical overview on the nomenclatural issues of the type species E. permagnificus are provided. Furthermore, this species is documented for the first time from Israel in association with Quercus calliprinos. The previously described North American species Exsudoporus frostii and E. floridanus are molecularly confirmed as representatives of Exsudoporus, and E. floridanus is epitypified. The eastern Asian species Leccinum rubrum is assigned here to Exsudoporus based on molecular evidence, and a new combination is proposed. Sequence data from the original material of the Japanese Boletus kermesinus were generated, and its conspecificity with L. rubrum is inferred as formerly presumed based on morphology. Four additional cryptic species from North and Central America previously misdetermined as either B. frostii or B. floridanus are phylogenetically placed but remain undescribed due to the paucity of available material. Boletus weberi (syn. B. pseudofrostii) and Xerocomus cf. mcrobbii cluster outside of Exsudoporus and are herein assigned to the recently described genus Amoenoboletus. Biogeographic distribution patterns are elucidated, and a dichotomous key to all known species of Exsudoporus worldwide is presented.

Evolutionary innovations through gain and loss of genes in the ectomycorrhizal Boletales.

We aimed to identify genomic traits of transitions to ectomycorrhizal ecology within the Boletales by comparing the genomes of 21 symbiotrophic species with their saprotrophic brown-rot relatives. Gene duplication rate is constant along the backbone of Boletales phylogeny with large loss events in several lineages, while gene family expansion sharply increased in the late Miocene, mostly in the Boletaceae. Ectomycorrhizal Boletales have a reduced set of plant cell-wall-degrading enzymes (PCWDEs) compared with their brown-rot relatives. However, the various lineages retain distinct sets of PCWDEs, suggesting that, over their evolutionary history, symbiotic Boletales have become functionally diverse. A smaller PCWDE repertoire was found in Sclerodermatineae. The gene repertoire of several lignocellulose oxidoreductases (e.g. laccases) is similar in brown-rot and ectomycorrhizal species, suggesting that symbiotic Boletales are capable of mild lignocellulose decomposition. Transposable element (TE) proliferation contributed to the higher evolutionary rate of genes encoding effector-like small secreted proteins, proteases, and lipases. On the other hand, we showed that the loss of secreted CAZymes was not related to TE activity but to DNA decay. This study provides novel insights on our understanding of the mechanisms influencing the evolutionary diversification of symbiotic boletes.

Evolutionary Morphogenesis of Sexual Fruiting Bodies in Basidiomycota: Toward a New Evo-Devo Synthesis.

The development of sexual fruiting bodies is one of the most complex morphogenetic processes in fungi. Mycologists have long been fascinated by the morphological and developmental diversity of fruiting bodies; however, evolutionary developmental biology of fungi still lags significantly behind that of animals or plants. Here, we summarize the current state of knowledge on fruiting bodies of mushroom-forming Basidiomycota, focusing on phylogenetic and developmental biology. Phylogenetic approaches have revealed a complex history of morphological transformations and convergence in fruiting body morphologies. Frequent transformations and convergence is characteristic of fruiting bodies in contrast to animals or plants, where main body plans are highly conserved. At the same time, insights into the genetic bases of fruiting body development have been achieved using forward and reverse genetic approaches in selected model systems. Phylogenetic and developmental studies of fruiting bodies have each yielded major advances, but they have produced largely disjunct bodies of knowledge. An integrative approach, combining phylogenetic, developmental, and functional biology, is needed to achieve a true fungal evolutionary developmental biology (evo-devo) synthesis for fungal fruiting bodies.

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi.

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE 'nests', or dense aggregations of TEs. Some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

Evolution of opposing regulatory interactions underlies the emergence of eukaryotic cell cycle checkpoints.

In eukaryotes the entry into mitosis is initiated by activation of cyclin-dependent kinases (CDKs), which in turn activate a large number of protein kinases to induce all mitotic processes. The general view is that kinases are active in mitosis and phosphatases turn them off in interphase. Kinases activate each other by cross- and self-phosphorylation, while phosphatases remove these phosphate groups to inactivate kinases. Crucial exceptions to this general rule are the interphase kinase Wee1 and the mitotic phosphatase Cdc25. Together they directly control CDK in an opposite way of the general rule of mitotic phosphorylation and interphase dephosphorylation. Here we investigate why this opposite system emerged and got fixed in almost all eukaryotes. Our results show that this reversed action of a kinase-phosphatase pair, Wee1 and Cdc25, on CDK is particularly suited to establish a stable G2 phase and to add checkpoints to the cell cycle. We show that all these regulators appeared together in LECA (Last Eukaryote Common Ancestor) and co-evolved in eukaryotes, suggesting that this twist in kinase-phosphatase regulation was a crucial step happening at the emergence of eukaryotes.

Mycology: Rediscovery of a lost model fungus highlights the origin of mycorrhizal symbioses.

Arbuscular mycorrhizae (AM) are the most frequent symbioses of land plants. By reisolating a long-lost fungus from nature, a new study cracks the genomics of an enigmatic fungal-cyanobacterial partnership and reestablishes a valuable model for understanding the AM symbiosis.

Gene family expansions and transcriptome signatures uncover fungal adaptations to wood decay.

Because they comprise some of the most efficient wood-decayers, Polyporales fungi impact carbon cycling in forest environment. Despite continuous discoveries on the enzymatic machinery involved in wood decomposition, the vision on their evolutionary adaptation to wood decay and genome diversity remains incomplete. We combined the genome sequence information from 50 Polyporales species, including 26 newly sequenced genomes and sought for genomic and functional adaptations to wood decay through the analysis of genome composition and transcriptome responses to different carbon sources. The genomes of Polyporales from different phylogenetic clades showed poor conservation in macrosynteny, indicative of genome rearrangements. We observed different gene family expansion/contraction histories for plant cell wall degrading enzymes in core polyporoids and phlebioids and captured expansions for genes involved in signalling and regulation in the lineages of white rotters. Furthermore, we identified conserved cupredoxins, thaumatin-like proteins and lytic polysaccharide monooxygenases with a yet uncharacterized appended module as new candidate players in wood decomposition. Given the current need for enzymatic toolkits dedicated to the transformation of renewable carbon sources, the observed genomic diversity among Polyporales strengthens the relevance of mining Polyporales biodiversity to understand the molecular mechanisms of wood decay.

Hallmarks of Basidiomycete Soft- and White-Rot in Wood-Decay -Omics Data of Two Armillaria Species.

Wood-decaying Basidiomycetes are among the most efficient degraders of plant cell walls, making them key players in forest ecosystems, global carbon cycle, and in bio-based industries. Recent insights from -omics data revealed a high functional diversity of wood-decay strategies, especially among the traditional white-rot and brown-rot dichotomy. We examined the mechanistic bases of wood-decay in the conifer-specialists Armillaria ostoyae and Armillaria cepistipes using transcriptomic and proteomic approaches. Armillaria spp. (Fungi, Basidiomycota) include devastating pathogens of temperate forests and saprotrophs that decay wood. They have been discussed as white-rot species, though their response to wood deviates from typical white-rotters. While we observed an upregulation of a diverse suite of plant cell wall degrading enzymes, unlike white-rotters, they possess and express an atypical wood-decay repertoire in which pectinases and expansins are enriched, whereas lignin-decaying enzymes (LDEs) are generally downregulated. This combination of wood decay genes resembles the soft-rot of Ascomycota and appears widespread among Basidiomycota that produce a superficial white rot-like decay. These observations are consistent with ancestral soft-rot decay machinery conserved across asco- and Basidiomycota, a gain of efficient lignin-degrading ability in white-rot fungi and repeated, complete, or partial losses of LDE encoding gene repertoires in brown- and secondarily soft-rot fungi.

Fruiting body form, not nutritional mode, is the major driver of diversification in mushroom-forming fungi.

With ∼36,000 described species, Agaricomycetes are among the most successful groups of Fungi. Agaricomycetes display great diversity in fruiting body forms and nutritional modes. Most have pileate-stipitate fruiting bodies (with a cap and stalk), but the group also contains crust-like resupinate fungi, polypores, coral fungi, and gasteroid forms (e.g., puffballs and stinkhorns). Some Agaricomycetes enter into ectomycorrhizal symbioses with plants, while others are decayers (saprotrophs) or pathogens. We constructed a megaphylogeny of 8,400 species and used it to test the following five hypotheses regarding the evolution of morphological and ecological traits in Agaricomycetes and their impact on diversification: 1) resupinate forms are plesiomorphic, 2) pileate-stipitate forms promote diversification, 3) the evolution of gasteroid forms is irreversible, 4) the ectomycorrhizal (ECM) symbiosis promotes diversification, and 5) the evolution of ECM symbiosis is irreversible. The ancestor of Agaricomycetes was a saprotroph with a resupinate fruiting body. There have been 462 transitions in the examined morphologies, including 123 origins of gasteroid forms. Reversals of gasteroid forms are highly unlikely but cannot be rejected. Pileate-stipitate forms are correlated with elevated diversification rates, suggesting that this morphological trait is a key to the success of Agaricomycetes. ECM symbioses have evolved 36 times in Agaricomycetes, with several transformations to parasitism. Across the entire 8,400-species phylogeny, diversification rates of ectomycorrhizal lineages are no greater than those of saprotrophic lineages. However, some ECM lineages have elevated diversification rates compared to their non-ECM sister clades, suggesting that the evolution of symbioses may act as a key innovation at local phylogenetic scales.

Mycena genomes resolve the evolution of fungal bioluminescence.

Mushroom-forming fungi in the order Agaricales represent an independent origin of bioluminescence in the tree of life; yet the diversity, evolutionary history, and timing of the origin of fungal luciferases remain elusive. We sequenced the genomes and transcriptomes of five bonnet mushroom species (Mycena spp.), a diverse lineage comprising the majority of bioluminescent fungi. Two species with haploid genome assemblies ∼150 Mb are among the largest in Agaricales, and we found that a variety of repeats between Mycena species were differentially mediated by DNA methylation. We show that bioluminescence evolved in the last common ancestor of mycenoid and the marasmioid clade of Agaricales and was maintained through at least 160 million years of evolution. Analyses of synteny across genomes of bioluminescent species resolved how the luciferase cluster was derived by duplication and translocation, frequently rearranged and lost in most Mycena species, but conserved in the Armillaria lineage. Luciferase cluster members were coexpressed across developmental stages, with the highest expression in fruiting body caps and stipes, suggesting fruiting-related adaptive functions. Our results contribute to understanding a de novo origin of bioluminescence and the corresponding gene cluster in a diverse group of enigmatic fungal species.

Transcriptome data reveal conserved patterns of fruiting body development and response to heat stress in the mushroom-forming fungus Flammulina filiformis.

Mushroom-forming fungi are complex multicellular organisms that form the basis of a large industry, yet, our understanding of the mechanisms of mushroom development and its responses to various stresses remains limited. The winter mushroom (Flammulina filiformis) is cultivated at a large commercial scale in East Asia and is a species with a preference for low temperatures. This study investigated fruiting body development in F. filiformis by comparing transcriptomes of 4 developmental stages, and compared the developmental genes to a 200-genome dataset to identify conserved genes involved in fruiting body development, and examined the response of heat sensitive and -resistant strains to heat stress. Our data revealed widely conserved genes involved in primordium development of F. filiformis, many of which originated before the emergence of the Agaricomycetes, indicating co-option for complex multicellularity during evolution. We also revealed several notable fruiting-specific genes, including the genes with conserved stipe-specific expression patterns and the others which related to sexual development, water absorption, basidium formation and sporulation, among others. Comparative analysis revealed that heat stress induced more genes in the heat resistant strain (M1) than in the heat sensitive one (XR). Of particular importance are the hsp70, hsp90 and fes1 genes, which may facilitate the adjustment to heat stress in the early stages of fruiting body development. These data highlighted novel genes involved in complex multicellular development in fungi and aid further studies on gene function and efforts to improve the productivity and heat tolerance in mushroom-forming fungi.